Shocks in Quark-Gluon Plasmas 
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Hydrodynamics is known to describe matter created in high energy heavy ion collisions well. Large 
deposition of energy by passing jets should create not only the sound waves, already discussed in 
literature, but also the shocks waves of finite amplitude. This paper is an introduction to relativistic 
shocks, which go through elementary energy and momentum continuity argument, to weak shocks 
treated in Navier-Stokes approximation, to out-of-equilibrium setting of AdS/CFT. While we have 
not yet found numerical solution to corresponding Einstein equations, we have found a variational 
approximation to the sum of their squares. Our general conclusion is that deviations from LS and 
NS hydrodynamical shock profiles are surprisingly small, even for strong shocks. We end with a list 
of open questions which the exact solution should be able to answer. 



PACS numbers: 



I. INTRODUCTION 



Shocks are rather well known hydrodynamical phe- 
nomenon described in the textbooks such as yy. Their 
discussion in the framework of ultrarelativistic collision 
of nuclei has been initiated by Landau and Belenky [2] 
who had applied it to matter compression at the initial 
time. Important point made in their paper was that very 
strong shocks (of very large amplitude) do not depend on 
the EOS of the matter before the shock. 

In 1970's, with the first experiments with ultrarela- 
tivistic nuclei at BEVALAC, it has been suggested to 
look for shock waves and Mach cones in excited nu- 
clear matter. Unfortunately hydrodynamics itself has 
not worked out in this application because the nuclear 
matter is not a good enough liquid: the nucleon mean 
free path is comparable to the nuclear sizes. Another 
applications of the shock theory to heavy ion collisions 
has been considerations of the deflagration and detona- 
tion shocks, propagating from the system's edge inward 
[Ij . Small velocity of such fronts were predicted to induce 
the so called "burning log" scenario for RHIC collisions. 
In practice however it was later found that a very long 
time needed for those shock propagation is not really 
available in real collisions, in which the matter outward 
gradients are large enough to generate strong radial flow 
overwhelming the deflagration. 

Unlike excited nuclear matter, Quark-Gluon Plasma 
(QGP) turned out to be a very good liquid, with the 
effective mean free path about an order of magnitude 
smaller than that for nucleons in nuclear matter. Soon 
after beginning of the RHIC era of experiments it became 
clear that already an ideal relativistic hydrodynamics de- 
scribes those very well [SHZj while the viscous corrections 
are indeed small ^8] |9] . Further recent confirmation came 
from the next-order applications, such as sound pertur- 
bations from initial state fluctuations. The calculated 
higher angular harmonics of the correlation functions are 
in good agreement with the RHIC/LHC data, see e.g. 

m- 



AdS/CFT correspondence has been used to explain the 
rapid onset of hydrodynamical regime. It started with 
the small viscosity-to-entropy ratio predicted in |llj 



rj/s = l/47r 



(1.1) 



for any strongly coupled theory with the holographic 
dual. Later those studies has been followed by studies 
of various out-of-equilibrium settings [I3HT5] . Basically 
all of them found rapid onset of viscous hydrodynam- 
ics, basically in the so called "infalling" time needed to 
approach the forming black hole horizon. 

Unfortunately, the out-of-equilibrium settings just 
mentioned are all time-dependent, and (except for the 
first one) require solution of the 2-1-1 dimensional dynam- 
ical Einstein equations, which is technically very chal- 
lenging. In this paper we propose to study the shock 
waves as an alternative out-of-equilibrium stationary set- 
ting. In this problem the Einstein equations are basically 
elliptic and thus allow for much easier treatment than the 
dynamical (hyperbolic) ones. In fact, we will be arguing 
below that one can use rather simplistic variational ap- 
proach and get sufficiently accurate approximation, even 
for strong shocks. (Weak shocks in AdS/CFT setting has 
already been discussed in Refs [TSJ [TH] and we will not 
study/review this case.) 

Another tool to be discussed is the so called resumed 
hydrodynam,ics 12 which suggests an approximate way 
to include higher gradients in sQGP. We will discuss 
the convergence of the corresponding series in gradients, 
as well as the corresponding Lublinsky-Shuryak (LS) re- 
sumption of those. 

Last topic in the introduction is phenomenological: it 
deals with some estimates for shock which appear due to 
jet quenching at LHC. We will argue that the energy de- 
position is large enough to create strong shocks at early 
time. The main point is the observation of very asym- 
metric events, in which the trigger jet has the transverse 
energy Et significantly larger than that of the associate 
jet Ea- Energy of the jet is deposited into QGP can be 



as large as 



AE fn Et - Ea ^ 100 GeV 



(1.2) 



(The first equality we write as approximate since the trig- 
ger also looses some energy, as well as picking up some 
from fluctuations and trigger bias effect, both of the scale 
10 GeV or so.) 

Let us now compare to the volume in which this en- 
ergy is deposited. Basic unit of length we use is the holo- 
graphic distance to the horizon l/{TrT) « .15 fm. The 
strong shock example we used above produces a width- 
at-half maximum of about 4 such units, making it about 
Sr « 0.6/to for T = AGeV = 1/(0.5 /m) . Energy of 
the jet is deposited into a cylinder of length L and radius 
given by the shock width. Its volume is 

V = Lnr^ - (10/m) * 3.14 * (0.36 fm^) « lO/m^ (1.3) 

If this energy AE is, say, 100 GeV, the deposited en- 
ergy density 6e ~ 10 GeV/frn^, comparable to the QGP 
pressure in the bulk 



PQGP~4T'^ = 4*(0.4Gey)/(0.5/m)= 
10 GeV/fm^ 



(1.4) 



The conclusions is obvious: such jet quenching event 
changes the stress tensor by some factor 2 or so. This 
should create what we called strong shocks, which re- 
main so at least for some time. What it means is that 
the Mach cones from them should be different from those 
considered in the literature [201 EI] : those were assumed 
to be sounds and move with the speed of sound. The dif- 
ference can be seen from the results of the nest section. 
The detailed study of such modified Mach cones can be 
given elsewhere. 



II. TWO CONTINUITY EQUATIONS 

For self-consistency of the paper we briefly review some 
well know material. 

Before we do so, let us introduce kinematical notations. 
The shock problem can be considered in two different 
frames. In the relativistic notations we will consistently 
use rapidity variables, in which Lorentz transformation 
between them is simply additive: we think it is pedagog- 
ically simpler than using Lorcntzian square roots. 

Frame A: the matter in front of the shock is considered 
to be at rest (zero rapidity). We denote by Y the rapidity 
of the matter behind the shock, and by y the rapidity of 
the shock itself. 

Frame B: The shock is at rest and the shock front is 
stationary. Shifting all rapidities by — y, one get inflowing 
matter with the rapidity Y~y and for out floating matter 
the rapidity is y. 

One difference of the QGP with the usual non- 
relativistic liquids is that the conservation of matter - 
of the number of atoms or at least some conserved vec- 
tor current - does not play an important role. Indeed the 



conservation of two vector currents available - that of the 
baryon number and strangeness - can well be ignored, 
as the values in question are very small. Therefore, the 
matter will be described by its stress tensor alone. The 
fundamental laws to be used are those of the energy and 
momentum conservation. 

There are no time derivatives and the vanishing di- 
vergences in hydro equations integrate to the continuity 
equations, of the T°i and T" : 

4:pfCosh{y — Y)sinh{y — Y) 

= ApiCosh[y)sinh{y) (2-1) 

3pfSinh{Y — y)^ + pfCosh{Y — yY 

= 3pisinh{y) + piCosh{y) (2-2) 

where we have used standard relativistic local expression 
for the stress tensor, with eij,Pij being the initial/final 
energy density and pressure on both sides of the shock. 

Let us start with the simplest conformal equation of 
state, the same QGP phase on both sides, which means 
that one can use e^ = 3pi. Now these two equations can 
be written as two expressions for the compression ratio 



El 
Pi 



cosh{y)sinh{y) 



cosh{y 
{3sinh{y)^ 4 



- Y)sinh{y 
cosh{y)^) 



Y) 



{3sinh{-Y + yY + cosh{-Y + y)2) 



(2.3) 



The second equality is in fact the equation for y, provided 
the rapidity jump Y of the matter is given. This equation 
has the following solution 



y{Y) =Y - (l/2)/n(e2^ + 1 - ^ e^^ + e^^ + i) (2.4) 

which is plotted, together with y{Y) — Y , in the upper 
part of FigjlJ 

The constant value of the shock rapidity at small com- 
pression (compression ratio close to 1, or weak shocks at 
F < 1) is 



y(Y^0) = {l/2)ln{2~V3) 
tanh{y{Y = 0)) = ^ 



(2.5) 



It is nothing else but the sound rapidity (velocity) . Thus 
two branches of the solutions are often called the super- 
sonic and subsonic matter flows. 

As the compression ratio grows, together with Y, to 
values of the order one Y ~ 0(1) one finds the so called 
strong shocks, to be discussed in what follows. Note that 
for initial state perturbations in heavy ion collisions the 
compressions are indeed generally of this category. When 
the compression ratio and rapidity jump gets large, those 
shocks can be called very strong. We will not discuss 
those: and just note that these two quantities grow very 
differently because of hyperbolic functions involved. For 
example large compression ratio of 10 corresponds to the 
rapidity jump which is still Y ^ I.. 




FIG. 1: (color online) (Upper plot) The rapidity of the shock 
y versus the matter rapidity Y is shown by upper (red) solid 
line, while the lower (green) dashed line is y — Y. (Lower plot) 
The compression ratio versus Y . 



This has consequences for the "Mach angle" as the 
cosOm — Vshock/c is approaching one. It is also important 
for the size of the "shock horizon" which is the distance 
the shock can travel till freezeout. 



III. WEAK RELATIVISTIC SHOCKS 

Weak, we repeat, means small rapidity jump through 
the shock front 



5^ = 2// - y^ < 1 



(3.1) 



In this case, as we will soon see, the gradients are small 
and this classic textbook problem is solved in the Navier- 
Stokes (NS) approximation (see e.g. [T] , chapter 87) . The 
principal steps in the relativistic case are the same as in 
non relativistic case, minor modifications are due to (i) 
the relativistic kinematics of the flow; (ii) the absence 
of the conserved matter current, and (iii) different EOS 
compared to the nonrelativistic gases usually considered. 



As we already stated, we ignore vector current and 
their conservations altogether, focussing on the stress 
tensor. Its relativistic dissipative part can be written 
as follows 



5T„, 



+e(A^^VpUp) 



-Ap^VpUp) 



(3.2) 



where the coefficients 77, f are the shear and the bulk vis- 
cosities, the projection operator onto the matter rest 
frame is 



^^U^f^ ^fiu — 5/iI^ 



(3.3) 



The 11 and 01 NS equations, for the simplest conformal 
EOS e = 3p, read 



pAuQUi + 7]dxUo ~ CIO 



(3.4) 
(3.5) 



where two constants in the r.h.s. can be inferred e.g. 
from flow (fluxes) far before the shock x — ^ —00 where 
the matter is homogeneous and the gradient terms are 
absent. Writing the functions as initial values plus modi- 
fications p(a;) = pi + Sp{x), y{x) = yi + Sy{x) and substi- 
tuting them to the two equations above, one can perform 
expansion in small terms up to the second order. Note 
that the viscosity term can be kept as constant as the 
gradient is already of the desired magnitude of small- 
ness. Since Sp only appears linearly, one can find it from 
one equation and substitute it into another, obtaining 
thus a closed differential equation for Sy{x) alone. It is 
quadratic in rapidity perturbation and can be rewritten 
in the following transparent form 



{Sy)i6y-Y) + {Ax) 

where, we remind, Y = yf - 
cient is (c^ = cosh{yi)) 



dSy 
dx 







(3.6) 



The value of the coefR- 



Ax = 



77 c,(cf-l)(4c| + 9) 
Up, -2cj + 6c| - 3 



(3.7) 



in which the initial rapidity can be approximated by that 
of the sound. 

Indeed, by construction, one root of the quadratic form 
is zero and the second root of the l.h.s. (other than y^) 
must be the jump to the final rapidity j//. This is the 
same generic equation as one gets (for pressure) of the 
non relativistic shock, and its solution is predictably the 
"Fermi step function" 



dy{x) 



Y 

1 + exp{—Yx/Ax) 



(3.8) 



The width of the shock for weak shocks is paramctrically 
large, 0{1/Y), which justifies the applicability of the gra- 
dient expansion and explains why it propagates with the 
speed of sound. 




FIG. 2: (color online) An example of a solution to the NS 
hydrodynamics, the rapidity (black dashed line) and pressure 
(blue solid line. Pressure is in units of its initial value: thus 
the curve starts near 1 on the left. 



IV. STRONG SHOCKS AND THE RESUMED 
HYDRODYNAMICS 

Let us now proceed to shocks with the rapidity jump 
Y = 0(1), for which there is no apparent small parameter 
and thus the NS equation cannot be expanded. Since 
pressure appears only linearly, one can easily manipulate 
01 and 11 equations into one single differential equation 
for rapidity (the second equality) 



P = 



-{4/3)Tjc{dy/dx)+M4s^,+l) 

_ -T]s{dy/dx) + ApiSjCj 
4sc 



where we use the short-hand notations 



(4.1) 



c = cosh{y{x)), s — sinh{y{x)) 

All quantities with the index i = initial are the cor- 
responding values before the shock, at x = — oo. The 
solution to the NS equation cannot be obtained analyti- 
cally, so we use a numerical solver. A particular example 
is shown in FigJ2j for the rapidity and the pressure. Its 
parameter were selected at random: the rapidity jumps 
up by a factor of 2 and the pressure down by about factor 
3. 

Now, how reliable is the NS solution, for such a strong 
shock? The implicit assumption is that all nonlocal 
(higher gradient) terms are small. We need to know 
(i) the individual coefficients - higher viscosities - of 
those terms; 

(ii) the values of the higher gradients; 
and (iii) the combined effect of their sum, or convergence 
of the series. 



For usual fluids such as water to air, the gradients are 
usually so small that we do not need higher viscosities. In 
practice their empirical values are not even known, and 
so it would be hard to even estimate the magnitude of 
those terms. Yet for sQGP we have the AdS/CFT cor- 
respondence, which provides in principle a complete set 
of such higher viscosities, and many of those has already 
been evaluated in literature. The sumnrary of these re- 
sults and their resummation was the subject of the "im- 
proved hydrodynamics" by Lublinsky and myself [I2] , LS 
for short. 

The linearized correlators in AdS/CFT can in principle 
provide the values of all kinetic coefficients, as the coef- 
ficient of certain powers of w and k (the frequency and 
the wave vector). Since there are several kinematically 
different channels excited by stress tensor, these calcula- 
tions also provide sufficient number of crosschecks. For 



example, ensuring that the first viscosity always is (1.1) 
anywhere viscosity appears. The "improved hydrody- 
namics" tries to combine all known coefficient into some 
"resumed" model functions. The main result was the 3- 
term Model 1 which is based on FADE approximation, 
which may be called LSI. It reproduces exactly eight 
first coeflacients and overall behavior of the correlators 
quite well. Furthermore, we found that the second and 
the third poles largely cancel each other, and for an es- 
timate of the effect of higher gradients we will use our 
simpler model LS2 in which the effective viscosity is 



1LS2 



Vo 



1 - V2.oky{2TrT)^ - zw77o,i/(27rr) 



(4.2) 



written in units in which all eta coefficients are dimen- 
sionless numbers. Furthermore, as there is no time de- 
pendence in stationary problems we now discuss, we will 
only need one simple coefficient 



'72,0 = 



(4.3) 



Due to its sign, the resumed factor reduces the effect 
of the NS term as k grows. This is indeed what was 
known for the lowest quasi normal modes (sounds): its 
imaginary part grows as k^ till some value, and then it 
stops growing. 

As the rest of this paper below will be discussing the 
AdS/CFT setting, it is convenient to us the holographic 
thermal horizon as the natural unit of length. Thus the 
horizon location in the holographic coordinate is 



Zh = I/ttT = 1 



(4.4) 



Using it and the holographic value of the first viscosity 
(1.1) one finds conveniently that the combinations which 



often appears is greatly simplified, for example the ratio 
in the r.h.s. of NS equation is 

rj rj A 1 
p s T ttT 



= 1 



(4.5) 



Of course, our problem has with two temperatures and 
pressures, pi and pf at both a; — > ±oo: but as the final 



0.03- 






0.02- 


A 

i/:^\ 




0.01- 


/ A 




0- 




h ^^ 






\ / 


0.01- 




\j/ 


0.02- 




w 


0.03- 







be the Fermi-like shape 



/ = 



0.03- 








0.02- 


fr 


\ 




0.01- 




■\\. 




0- 






y ^ 


0.01- 






\ 7 


0.02- 






w 


0.03- 









-2 



FIG. 3: (color online) (a) Gradient expansion for the Fermi 
step with the width w = 0.75. By the (red) dashed, (blue) 
dash-dotted, (black) solid and (brown) dotted lines we show 
the results including up to 2,4,6 and 12 gradients. The lower 
plot shows the LS modification factor 1/(1 — y^' /8y) applied 
to the same shape, by black solid line. (Red) dash-dotted and 
(brown) dotted line have up to 8 and 12 derivatives, respec- 
tively. 



parameters do not appear explicitly in the equations in 
the form used, we use the initial temperature to set the 
units. 

The first issue to be discussed is the magnitude of the 
higher gradient terms and convergence of the derivative 
expansion. Using for simplicity our LS2 model we write 



1 + 978 + (978)^ + 



(4.6) 



1 + exp{~x/'w) 



(4.7) 



1 - (92/8 
and apply it to the shock profile, taken for illustration to 



Weak shocks have large width w >> 1 and thus small 
gradients, so that the convergence is parametrically jus- 
tified. Proceeding to strong shocks, such as our example 
displayed in Figj2j one finds the width of the order one: 
thus no small parameter is available. Few terms in the 
expansion for this example are shown in Fig[3J The first 
thing to notice is the fact that all corrections are below 
1 percent, more than order of magnitude below the level 
which the numerical smallness 1/8 of the expansion sug- 
gest. 

To investigate the issue further, one may reduce the 
width and check how this expansion in gradients be- 
haves. We observed the first signs of trouble at the width 
w — 0.75. As seen from the upper Fig|3J the first terms 
indicate a reasonably small correction, with the mag- 
nitude of 2%: but going further (see the brown curve 
with 12 derivatives) one finds widely oscillating correc- 
tions of increasing amplitude. This is a behavior typical 
for asymptotic series, which converge only till the certain 
term. It is usually assumed in such cases that the last 
"good" term shows the best possible approximation: we 
take the black curve including up to 4th correction (8-th 
derivatives) to by the closest to the truth. 

The next issue is if the "resumed" version in which the 
first correction (two derivatives) is put into the denomi- 
nator. In other words, the question is if one can naively 
sum the geometric series, ignoring the operator nature of 
the derivative. (Of course, the exact resumption into the 
denominator is provided by the inverse of the differential 
operator in question.) The result is shown in the lower 
FiglS] by the black curve, which is indeed quite close in 
shape to the sum of the first 4 corrections (the red curve) . 

We have also applied the same operations as for the 
Fermi step functions above to the numerical solution of 
the NS equation (4.1) shown in Fig[2J The results are 
very much the same as above for the Fermi step func- 
tion. The upper figure shows what happens with a series 
in subsequent gradient expansion: it seems well conver- 
gent first, but the 6-th order term (12 derivatives) gets 
oscillating. The lower curve shows the LS resumed fac- 
tor, with the second derivative in the denominator. Note 
that the whole scale of effect is very small, about 0.003: 
we thus conclude that for practical purposes our strong 
shock example shows only negligible effect of the higher 
gradients. The NS approximation seem to be justified for 
this example, which was one of our main results. 

The conclusion from that is: the LS resumed hydrody- 
namics predicts, that the sum of the high gradient terms 
is much smaller than individual ones and basically van- 
ish. The shock is thus one more example of unexpectedly 
early onset of the NS behavior. 

It is still somewhat surprising how small is the LS re- 
sumption correction, even for a relatively strong shock 
example under consideration. Part of the explanation is 
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FIG. 4: (color online) The upper plot shows several terms of 
the gradient expansion for the NS solution shown in Fig[2] 
By the (red) dashed, (blue) dash-dotted, (black) solid and 
(brown) dotted lines we show the results including up to 2,4,6 
and 12 gradients. The lower plot shows the LS modification 
factor 1/(1 — y" /8y) applied to the same numerical solution 
of the NS equation: its magnitude and shape should be com- 
pared to the solid curve in the upper figure. 



perhaps the observation that the second derivative van- 
ishes at the inflection point, which is located right at the 
center of the shock solution. Such explanation would not 
be however applicable to two more recent examples of 
rapid onset of the NS regime discussed in Refs. [UJ [TS] . 



V. SHOCKS IN THE ADS/CFT 



We use coordinates v,xi — x,X2,X3,r and write the 
nonzero component of the metric as 



2 (■ 2 I 2 2 

ffii = -r fc +r s ; 

gi2 = g2i = -r^fcs + r^cs + A{xi,r)- 

922 = -r^fs^+r^c^ + B{xi,r); 

9i5 = 951 = c; 

925 = 952 = s; 

944 = 333 = r'^ 



(5.1) 



where 



/ = l-hix,r)yr^; 

c — cosh[y{x)]; s = sinh[y{x) 



(5.2) 



Note that we do look for a static (w-independent) solu- 
tion, depending on only 2 spatial coordinates r, x. The 
metric contains 4 unknown functions, y, h, A, B. 



If the two functions h and y are just constants, the 
Einstein equations (with appropriate cosmological con- 
stant) are satisfied without extra correction, namely with 
A,B = Q. Indeed, in this case this metric is nothing else 
but the black brane moving with the rapidity y. The 
shock is an interpolation between two such solutions, 
with yi , hi different from yf , hf in a finite- width region, 
the shock. Needless to say, the values are not arbitrary 
and must be related by the continuity of the energy and 
momentum flux, as we discussed at the beginning of the 
paper. 



The corresponding Einstein equations can be easily de- 
rived (e.g. by the "tensor package" of Maple), but they 
look very long and discouraging. Here is the simplest of 
them, the Ricci scalar 



R = - {-96 A"^ h^c^h"r'^- 96s Ah^c^h"r^B-mr^ 
+16scA^r*A" - Ur^B + 4r^A^ - 4AA"r^ 
+16c^h'^BA"A + 32c^Br^h^{h') 
+16c^B'^h^{h')r - UcBsr^A' 
-UB'r'^ + 2Ac^B^h^h"r^ 
+2c%''B^B" - 8AA"r^c'B + Acr'^A' 
+16r'^sA'^yc^ - Ar^c^B^sy 
+48sA^cA'r^ + 52cAB'r^s + Sscr^h'^A' 
+8c'^r^sAA + Wc'^rh^AA' - 'ic^rBh'^B' 
-l6c^rAA'h'^ - 52AA'r^c^B - 6AA^h!^c^{h')r 
+6AA^h^c^{h')r + 12sAcBr'^ - SsAh^r'^c 
-l6sAe'Bh'^ ~ 2r^B" + 3r^f - 3r^A'^ 
+16A^^ - 16A^c^ - l6sc''h^A^A" 
-AsAr^y - 6c^BA''^r'^ - lOc^Bfr'^ 
+20c^Br^y + 2c^Bfr^ + 16A^c^sB 
+2Ac''B^h^{h'f - Ac^h^BA'^ - c^h'^BB'^ 
+Ac'h^BA'^ + 2sr^A' - r^c'^h^B'^ 
+Ar'^c^h*A'^ - 16rA^Ac^ + 16r^d^BA 
+16A\''c^ + r^c^BB'^ - Ac^h^'^A'^ 
+AAW^A''' + 16A\''c'f 
-AA^r'^c^f - ASrc^A^A' + A?>rc^A^A' - 
Ar'^c'^A^A'^ - 12r2c6A2y2 + iGcVA^i 
+Se'rB'^A + Ar^e'B^y - Sr^c^B^B' 
-Sr^c^B^y^ _ Sr^c'^B^A'^ 
-~96r^c'^BA^ + Sr^c^BB' - r^c^B'^f 
+16r^d^AB - mA^r^B'c^ - Ac^r*yA 
+Ac^r^yA - 2c^r'^A'B 
+12c^B'^yr^ - 32c^r^A'A - SA^yr^c^ 
-2r^c^yB' + 2r^cyB' + ?,2r^c^AA' 
-32sc^ABh^{h')B' - 32r^sAc^h^{h')B' 
-32r^c^BsA'h^{h') - 2r^c^BAsyB' 
-288sAh^c^{h'fr^B + l6c\^Bh^{h')B' 
+6Ac'^r'^AA'h^h') - 6Ac\^AA'h^{h'^ 
+Ac^Bsyr'^A' - lAAsAh^c{h 



2„/u'\2 4 



r6AsA'^c^A'h^ih') - 6Ac'^BAA'h^{h') 
-6Asc'^A^A'h^{h') + Asc^h^BA'B' 
lAAsAe'B^h^ih'f + 6Ac^BAA'h^{h') 
-l6&>B^h^s(h')A' + Ar'^c^BsyA' 
-2r*c^BsA'B' + Sr'^c^BAA'y - 2Ac^rsAAB 
+2r^c''BAA'B' + 2r^c''B^syA' 
+Ar^c^BsA'A + 2r^&'BsyB 
-Ar'^e'BAA'y + Ar*sAc*yB' 
-6r'^sAc^yB' - '^r'^sAc^yA 
-Sr^sA^c^A'y - 16 sr^h^c{h') A' 
+12se'r^AfB + AAsc^r^ABB' 
+Asc^r^AA'B + Ssc^r^ABA'^ 
+8se'r'^AyA - Asc^A^r'^A'B' 
+56sc^A'^rA'B + Asr^h'^c^A'B' - 16sc'^r^BAy 

(5.3) 



28c'^BsAyr^ - Ar^csy + 2r^c^h'^B" 
16r^A^A'^ - 2r^c''B^B" - IGr^c^A^A" 
96A'^h^c^h"B + Sr^h^h" - Ar^B"c^B 
ASsAh'^ch^r'^ - ASsAj'B^h^h" - %c^h'^BB"sA 



Vic^h^J^B" + 8r^e-BB"sA - ^sC'h^BAr 



II Ji 



-Wr^c^ByA - I6c^ r"^ h^ A" A + 96A^h^c'^h"r 
+16c^r^yA - 16c^r'^yA + 8c^A^B"r^ 
+16c''r^h^A"A + 22A^ A^r^ - Sc'^h^A^B" 
+8c^B^h^h" - Wsr^c^BA' - l6sr^c^AB' 
+Asr^yA'c^ + Sr^Ac^A'A + Ar^Ac^yB 
-Sr^e'AA'y + Ar'^c^sA'A + 2r^c^syB 
-Sr^c^AA'i - Ar'^c'^AyB + Ar^c'^AA'B' 
-2r*c^AsB'^ - Src^A^B's + Arc^A^B'B 
+Ar'^c^A^yB' + 96r'^sAc!^B'^ - 8r^c*BsB 
-Ar'^c^ByA - 2r^c^BA'B + 2Ar^c*BAA' 
-16rc^B'^AA' + 8c*rA^Bs - Ac^rABB 
-Usc^r^B'^A' + Sr^AA'cy - 2sr^cA'B' 
+Ar'^c^ByA - S2r^sAAc^ - 16A^rc'^ys 
+8A^rc^yB + 22Ar^c^BAs - 2r^c^ByB' 
-AA^r^c^yB' + 8r^c^h!^{h')B' + Ac^r^AA'y 



y 



2A^A'r^B' -288A^h^c^{k 



2(,2„2^i,/\2„2 



+288^2/^2 c4(/i')V2 + 2c^Byr^B' ~ Schh^AA'^ 
+192A^h'^c^h')'^s - 288A^h^c'^{h')'^B + 32c^A'^h^h')B' 
+8c'^sh^AA'^ + 2c^sh'^AB'^ + 288c^BA^h^{h'f 
+8c^B^h^{h')B' - 8c^h'^AA'B' - l92sc^A^h^{h'f 
+8c''h^AA'B' - 32A2/i3c4(/i')B' + 8sAcA'\^ 
+20sAc^y V - AOsAc^r^y - 8sAcfr^ 
+72c''B^h^{h'fr^ + 12c^Bh''{h'fr^ 



3/1,' 



-6AsAe-Bh^{h 



\r + 8sr 



,6d" 



B"cA + l6r'^c^ByA 



+l92sAcr^ + l2cBr''y + Ac^Bh'^r^ 
-16scr^A' - 20c^r^AB - 8c^r^sB 
-20c^r^BB' + 6Qc^A^B'r^ - Ar^c^h^B' 
+96A^r'^c^B - 6AsAr^h^{h')c + 8sc^h^BA'r 
+8c^sh'^AB'r + 6AA^h^e'h"s - Ase'h^B^A" 
-l6c''h^BA"A - 7c2r6y2 + ^6^2^,2 

+2Ar'^h^{h'f - 8cr^A - AA^c^B^ 
-56r^c''B^ - 2Ar''c'B^ - 22AA^r^c^ 
+16r2c5^2^' - 8r^c^A^B" + IGr^c^yl^A" 
^Ar'e'B'^k + 8e'r^kB + Ac^h^Br^B" 
+16r^/i^(/i') + 16r^cy + 16cV'^i 
+8c^r^B' + \6A^c\^ - WA^c'h^ 
-12AA'r^ - 96c^Br^ - 36c^B^r'^ + Ac^B^h'^ 
+2Ac^Bh^h"r^ - GAsc'^A^h^h" - 16r^c^sA'A 
-~l6sc^r^AA'B + 96c^BA^h^h" - Wse'r^A^y 
-Ascr^h^A" + l6sc^A\'^A"B - 8c\sh^Ar''B" 
-8r'^c^syB + l6s&'h'^A^A" - Ar^c^B^AA") 

1 
2r2(r2-2csA + c2B)3 



(5.4) 
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FIG. 5: (color online) The scalar square of the modified Ein- 
stein tensor ( |5.7| l (upper plot) and the simple sum of squares 
of all components (lower plot) the as a function of the coor- 
dinate X (in units of l/nTi), for r = 1. The dashed (blue) 
curve is for the case without the B correction (5.81, the solid 
(black) line includes it . 



where the prime stands for the derivative over r and the 
dot means the derivative over xi. Substituting again 
the moving brane solution one gets constant (coordinate- 
independent) value R = 20, related to the AdS cosmo- 
logical constant. 

One can start by taking drastic simplifications: putting 
extra functions to zero A,B — and h{xi,r) = h{xi) 
to be independent on r. This results in the following 
relatively simple expression 



R-20 



1 
2t^ 
+4cosh{y)sinh{y)y 



— 16rcosh(y)y + y {7cosh{y) —3) 



(5.5) 



depending only on the rapidity profile across the shock. 
Let us see what this combination is for traditional Fermi- 
step solution in weak shock Example 1. Since in this case 
the width of the shock is large, the square of the deriva- 



tive and the second derivative are negligibly small, while 
only the first term is important. With this lesson in mind, 
one can return to a case with nonzero extra functions and 
use linearized equations, ignoring higher derivatives in x. 
This procedure leads to hydrodynamics, as one can solve 
all the equations for extra function perturbatively, see 
[T7] for general review of the method and Refs [THl |T5] 
for the particular case of weak shocks. 

Let us turn to strong shocks, such as our numerical 
example, in which case the width of the shock is 0(1) 
and there is no small parameter in the problem. The 
Einstein equations in full nonlinear form take pages and 
are way too complicated to be presented here or solved 
directly. Instead we propose to use a variational method. 

Unfortunately the Einstein-Hilbert action R is not 
bounded from below and cannot be used for variational 
studies. The so called conformal gravity, with a squared 
Weyl tensor in the Lagrangian, should work [5?. What 
we propose to do is to use the covariantly squared (mod- 
ified) Einstein tensor 



E' 






E„ 



Emn + 65,, 



(5.6) 



which combines all the Einstein equations (in the 
AdS/CFT setting) into one (covariant scalar) combina- 
tion. For a check, one can of course look at all individual 
components of the modified Einstein tensor. We used for 
this purpose a (non-covariant!) sum of squares 



ES — y ^{Efnn) 



(5.7) 



This combination is of course sign-definite, and since for 
a solution all components, and thus the sum of squares, 
should vanish, it can be used to monitor the "variational 
progress" for all of the components. 

The equations are basically elliptic, and as such they 
only need the field values on the boundary of the region 
to be appropriately interpolated (solved) inside it. In this 
problem all corrections vanish both at a; — >■ ±cxd, as well 
as at large r. The only tricky issue is near the nontrivial 
boundary of the black hole horizon. We use as the initial 
input the rapidity y{x) and pressure p{x) = h'^{x) from 
the numerical NS solution, which fix the gyy, but allow 
nonzero modification functions in other components such 
as g^x ■ The inter mediate step we will not dwell on was a 
parameterization of the solution in simple enough form, 
so that the evaluation of the huge expressions involved 
be possible. 

A simple way to proceed is to use variationally, to use 
certain ansatz (assumed trial function) and substitute it 
into all the equations and/or E'^ see how close/far are 
the results from the desired zero values. Evaluation of 
a not-too-complex trial function is performed by Maple 
in seconds, in spite of horrendously complex expressions 
involved. In FigJ5ja) one can see a comparison of the 
E'^{x,r) for the NS profile only (the top blue curves) 
with the results including this B (the lower black curves) . 
Note first, that even the dashed curve which use only 
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the input rapidity and pressure hydrodynamical solution 
(shown in FigM is aheady not too bad. Indeed, the 
components of the Einstein tensor are expected to be 
0(1), as the problem has no parameters, and the number 
of nonzero ones summed up is 11 (9 from the t,x^,x^ 
block plus E33 — i?44, other 6 are zero due to parity of 
x^ and X*). Thus one expects something of the order 
0(10) and gets something 20 times less at its peak. 

We have used only one of the correction functions B, 
and after some number of trials we came up with the 
flowing ansatz for it 

B{x, r) = -0.052r(l - O.Sx)exp[~.l{x + 0.3)^] (5.8) 

It is clear that the mismatch is reduced by another order 
of magnitude, in the whole region of x, r in question. 
(We don't show it as a function of r, but the relative 
error is constant in r). One can view this mismatch as 
the amount of external matter (or, more exactly, stress) 
which is needed to make our approximate solution of the 
Einstein equation exact. 

In the next Fig[5](b) we plot the simple sum of squares 
of all components of the (modified) Einstein tensor. 
While we have not minimized this quantity, it also shows 
a significant reduction of the mismatch. We take it as the 
indication that there are no significant cancellations be- 
tween components in the scalar and we are in fact close 
to the solution, to which one can get even closer with 
more sophisticated functions used. The standard proce- 
dure would be to discretize the functions, by introducing 
a grid in Xi,r and use well known relaxation methods 
(solving for zero) at each point [23] : in this case one can 
reach arbitrary high precision if needed. 

Of course, in the calculation we also were looking at 
the individual components of the modified Einstein ten- 
sor. Few are shown at Fig. |6j the reader can see that 
while in some cases ( e.g. the 15,33 components) the 
improvement is obvious, other components do not show 
that. The reader then is invited to look at the scale of all 
the graphs and notice that "bad" ones are much smaller 
than "good" ones naturally. Another way to see that is 
to note that the difference between the dashed and solid 
curves is due to correction B which was found to be only 
0.05 in magnitude, while the metric components are gen- 
erally not small 0(1). 

The meaning of the correction B is seen from the fact 
that it modifies the length elements along the coordinate 



dl 



rV-fs^ + c2 + B{x, r)/r^dx (5.9) 



Since B/r^ ^ 1/r one finds that far from the black hole 
the correction is unimportant and thus at the bound- 
ary this correction disappears. Near the horizon, since 
the first term is small, / « the second is dominant 
and 0(1). As the correction B is negative, it shrinks 
a bit the distance across the shock. The magnitude of 
B found leads to our main conclusion, that B makes a 
shock few percent sharper near horizon, as it is at the 
AdS boundary. 



(Additional comment about a horizon. In general, 
as new solution is found one has to calculate the null 
geodesies in it and find their bifurcation. Note however 
that for "radial" ones, with dx=0, the sign of the dr/dv 
is defined by the metric components 511, (?i5, which do 
not include the modification function B. So the line of 
horizon is still given by the (lorentz transformed) h'^{x) 
line.) 



VI. 



SUMMARY AND DISCUSSION 



Hydrodynamics is an approximate effective theory, its 
local approximation relies on smallness of the gradients. 
Another way to explain the assumption: the spatial scale 
of interest is much larger than the "micro scale" defined 
by the viscosity. (If kinetic theory be applicable, the lat- 
ter would be the mean free path of quasiparticles: but 
in our problem we have neither quasiparticles nor kinetic 
description.) Very weak shocks are basically sounds, the 
gradients are parametrically small and for them hydro- 
dynamical assumptions are justified. 

Strong shocks in general have no reasons to be treated 
in the NS approximation: there are no small parameter 
available and the higher gradient terms are in general not 
small. And yet we found, from two different (although 
related) sources, that all corrections to NS are of the 
magnitude of few percents or less. 

While for most applications one usually don't even 
know the coefficients of the higher gradients, in the case 
of conformal plasmas the AdS/CFT correspondence tells 
us their values, and several of them has been calculated. 
Lublinsky and myself [12] (LS) had proposed to resum 
those into some "universal hydro" form. The alternating 
signs of the corrections basically lead to strong cancel- 
lation of those higher order terms, even when they are 
not small individually. Explicit calculation of the higher 
derivative terms show that the series are asymptotic. The 
LS resumed correction indeed provides results which are 
very consistent with the convergent part of the gradient 
expansion. They are also small, indicating that solutions 
to "universal LS hydrodynamics" are remarkably close to 
the NS ones, at least in the case of shocks. 

The second approach used is the AdS/CFT. While 
looking for shock-related stationary solutions to the Ein- 
stein equations we made the first steps, by using some 
original form of the variational approach. We have been 
able to reduce that scalar square of modified Einstein ten- 
sor to a fraction of a percent, in the whole range of vari- 
ables, which perhaps indicate we are numerically close to 
a true solution. Unlike earlier studies [IH] , there is no 
expansion in gradients or linearization of the equations 
in our approach: all terms are kept and all of them are 
individually 0(1). 

Our variational solution finds that one needs only a 
few percent correction to such observables as the shock 
width and profile, as compared to the one obtained from 
the NS equations. Our results agree in spirit with finding 
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FIG. 6: (color online) The same notations as for the previous plot, for individual components of the modified Einstein tensor, 
from top to bottom £11, £22, £15, Ezz 



from two more recent examples [HllTSj of the onset of the 
LS/NS hydrodynamics. Those two examples are however 
time-dependent collision problems, which are much more 
complicated technically than the one described in this 
work. 

Needless to say, a lot of work needs to be done. The 
variational solution should be extended to a better accu- 
racy. The issue may be pursued more into the realm of 
even stronger shocks. Last but not least, one now may ar- 
gue that production of such shocks by strongly quenched 



jets at RHIC/ LHC heavy ion collisions opens a possi- 
bility to test some of these predictions experimentally. 
The first steps perhaps would be observation of the cor- 
responding Mach cone angles in jet-hadron correlations 
and derivation of the shock wave velocities. 
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